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Computer Programs 


One major acoustic computer program and one performance preprocessor were developed under 
the counterrotation portion of this Contract described in Volume I. 

CRPFAN (CounterRotation PropFan Noise) predicts 

a. Steady loading and thickness noise of Rotor 1. 

b. Steady loading and thickness noise of Rotor 2. 

c. Unsteady loading noise of Rotor 2 resulting from the interaction of the blades with wakes 
and vortices shed from Rotor 1. 

d. Additional BPF harmonic noise resulting from the interaction of the rotors with non- 
uniform inflow generated by the installation environment (see Volume I, Section 2.2). 
The installation effects considered are: angle-of-attack operation; the presence of a 
pylon; the presence of a fuselage, and the distortion behind a wing lifting line. This pre- 
diction is added to the results from (a) and (b) above. 

The operation of the code is sketched in Figure 1 . Depending on requirements, the program can 
be used to predict 

Items (a) - (b); 

(a) -(c); 

(a) -(d) 
or (d) alone. 

Three input files are required to run the program: 

Unit 10 Contains the geometric parameter for Rotor 1 together with universal 
operating conditions (for example temperature, pressure, pitch change 
axis spacing) and such control switches as are required. 

Unit 11 Contains unique geometric parameters for Rotor 2. 

Unit 12 Contains specific installation effects input. 

All input parameters are documented in the program. 

Two files are output: 

Unit 16 Contains the output from the “isolated” program. 

Unit 14 Contains the output from the installation effects program, including the 

combined “isolated + installation effects” results if desired. 

The output parameters are documented in the program. 

CRPCLD (CounterRotation Propeller Coefficients of Lift and Drag) computes blade section lift 
and drag coefficients as functions of angle-of-attack given counterrotation propeller performance 
map information as input. This information is used as installation effects input to CRPFAN. Three 
input files are required to run the program: 

Unit 10 Contains geometric parameters for Rotor 1. 

Unit 1 1 Contains geometric parameters for Rotor 2. 

Unit 12 Contains performance parameters. 
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All input parameters are documented in the program. 

Two files are output: 

Unit 6 Contains the full output in printout form. 

Unit 16 contains the calculated coefficients in “NAMELIST” form, suitable for 
inclusion into the input file for CRPFAN. 

The output is documented in the program. 

This Volume Contains: 

• Acoustic Program 

— CRPFAN Source Code 

— Sample CRPFAN input files: 

Rotor 1 geometry + operating conditions 

Rotor 2 geometry 
Installation Effects 

— Sample CRPFAN output files: 

Isolated 

Installation Effects 

• Performance Program 

— CRPCLD Source Code 
— Sample CRPCLD input files: 

Rotor 1 geometry 

Rotor 2 geometry 
Performance parameters 

— Sample CRPCLD output files 
Printout 

Namelist 


£ 
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Figure 1. Methodology - Counterrotation Prediction Model. 


3 










CRPFA N Pr ogram 
Description 



COUNTER-ROTATING PROPFAN NOISE PREDICTION PROGRAM 

CRPFAN 

The Main Program, CRPFAN, controls program execution. In addition, aerodynamic 
parameters and uninstalled tone noise are calculated for the two rotors, and printed as follows. 

Following the output of input quantities for both rotors (with input YMC & ZMC converted to 
Face Alignment and Mid-Chord Alignment) Acoustic predictions are output versus observer and 
emission angles for each rotor. 

» Tone predictions are identified by frequency and the No. ‘M’ as follows: 

M = 1 -> NHM: SL+T harmonics of Rotor 1 BPF 

NHM+1 -> 2*NHM: SL+T harmonics of Rotor 2 BPF 

2*NHM+1 -> NWHM*NHM+2*NHM: Rotor 1/Rotor 2 UL interaction tones, 

identified as nnBPFl+mBPF2 where: 
(M/NHM) = {(nn+1) + (m/NHM) >. 

For example: NHM=10 (recommended); NWHM=15: 

M = 1 -> 10; 1st 10 harmonics of Rotor 1 BPF 
1 1 -> 20; 1st 10 harmonics of Rotor 2 BPF 
21 -> 170; Rotor 1/Rotor 2 interaction tones: 

21= 1F+1A 
22= 1F+2A 

23= 1F+3A 

121 = 11F+1A ...etc. 

Finally, these tones are combined into 1/3 octaves for output 

The following Subroutines are called from CRPFAN: 

READIN - Reads input files 
GAMCAL - Calculates ratio of specific heats 
PITCH 1 - Calculates geometric parameters for Rotor 1 

LSPFIT - General interpolation and integration routine 
PITCH2 - - Calculates geometric parameters for Rotor 2 
RWVM - Rotor Wake/Vortex Model 
r SKCAL - Calculates Sears Function (Unsteady Lift Response) 

BJFN - Calculates Bessel (J) Functions 
> AKVCAL - Calculates chordwise non-compactness factor (thickness) 

AKLCAL - Calculates chordwise non-compactness factor (loading) 

TOBSPL - Combines tones into 1/3 octave bands 
ATMOSP - Applies correction for atmospheric absorption 
SPLOUT - Prints 1/3 octave output 

CRPIE - Counter Rotating Propfan Installation Effects model 
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SUBROUTINE READIN (Page 42) 

This Subroutine reads the input Namelist files containing geometric information for each rotor 
and program control switches (Namelist ‘INPUT’); and, (if the installed calculation is to be 
performed) .installation effects input (Namelist ‘INSTAL’). It rearranges the input geometry for 
input to the installation effects routine (CRPIE) and prints this input 

The following Subroutines are called from READIN; 

GAMCAL- Calculates ratio of specific heats 

LSPF1T - General interpolation and integration routine 


CRPFAN INPUT VARIABLES 


Namelist ‘INPUT’ 

Program Control 

ISOLAT - Controls Uninstalled CVR Prediction. 

=1, Compute “Isolated” Noise, Comprising; 

SL+T Noise of Rotor 1; 

SL+T Noise of Rotor 2; 

UL Noise of Rotor 2 Resulting from Interaction with the Wakes & Vortices 
shed from Rotor 1. 

=0, Do not perform Calculation. 

INST - Controls Installed C/R Prediction. 

=1, Compute “Installed” Noise resulting from: 

Angle-of- Attack effects; 

The presence of a Pylon; 

The presence of a Cylindrical Fuselage; 

The Presence of a Wing Lifting Line. 

For Rotor 1 & Rotor 2 harmonics of BPF. 

=0, Do not perform Calculation. 


NOTE: 

If INST=1 & ISOLAT=l the results from the two predictions are combined. 
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Definition of “Test” Conditions 
NCASE - No. of Cases 

DT - Dt; Tip Diameter, ft. 

HTR - Hub-to-Tip Radius Ratio. 

NB - B; No. of Rotor Blades. 

SHP - Shaft HorsePower. 

RPM - Rotational Speed, rpm. 

DELBP(*) - Pitch Angle Change from Design, Delta(Betap) deg. 

BETA34(*) - Reference Pitch Angle, Beta3/4 deg. 

Z34(*) - Actual Radius Ratio of Beta3/4. 

VO - Flight Speed, Vo ft/sec. 

TAMB - Ambient Temperature To deg. F. 

PAMB - Ambient Pressure Po psi. 

NHM - No. of BPF Harmonics to be computed. 

NSL - Sideline Indicator 

=1 Sideline prediction 
=0 Constant radius arc prediction. 

NTH - No. of Observer Angles, Theta-v. 

THAD(+) - Table (17 max.) of Observer angles, Theta-v deg. 

DIST - Sideline Distance ro ft. 

XPCA - Axial Distance between Pitch Change Axes, normalized by Rotor 1 Tip Diameter. 

NOTES: 

*If BETA34 is included, the value of DELBP in the NAMELIST input is overwritten. 

If Z34 is not present, a default value of 0.75 is assumed. 

+A default table is present in the program. 

Definition of Blading 

NZ - No. of Spanwise Stations 

Z(+) - Radius Ratio; z=r/rt 

ZMC(+) - Axial Location of Blade Mid-Chord Line relative to 
Pitch-Change Axis, normalized by Tip Diameter, Dt. 
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YMC(+) 

CHORD(+) 

TMOC(+) 

BETAP(+) 

THETA(+) 


- Tangential Location of Blade Mid-Chord Line relative 
to Pitch-Change Axis, normalized by Tip Diameter Dt 

- Blade Chord normalized by Tip Diameter, b/Dt. 

- Maximum Blade Section Thickness, normalized by Chord, tm/b. 

- Design or Reference Blade Section Pitch Angle, Betap, deg. 

- Blade Section Swirl Coefficient, Theta=(Cu/U)exit. 


NOTES: 

+NZ required, maximum of 10. 


Steady Loading and Thickness Options 

ITDO - Thickness Distribution Option Indicator: 

=0, Point-Impulsive. 

=1, Rectangular. 

=2, Parabolic. 

=3, Triangular. 

=4, Generalized Trapezoidal. 

=5, Generalized Parabolic. 

ALE++ - Leading Edge Thickness normalized by Max. Thickness, ale/tm. 

BTE++ - Trailing Edge Thickness normalized by Max. Thickness, bte/tm. 

XTM++ - Position of Max. Thickness relative to Mid-Chord, normalized by Chord. 

DXTM - Width of Max. Thickness Region normalized by Chord. (ITDO=4 only). 

ILDO - Loading Distribution Option Indicator: 

=0, Point-Impulsive. 

=1, Rectangular. 

=2, Linear. 

=3, Generalized Parabolic. 

XLM - Position of Max. Loading relative to Mid-Chord, normalized by Chord 

(ILDO=3 only). 


NOTES: 

++For ITDO=4 and rrDO=5 only. 

General Input 

IDS - Doppler Shift Indicator 
=0, No Doppler Shift. 
=1, With Doppler Shift. 
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IPTO - 1/3 - Octave Band Printout Indicator 
=0, No 1/3 - Octave Printout. 

=1, With 1/3 - Octave Printout. 

ISECT - Interpolated Spanwise Sections Indicator 

=0, Calculations performed using input stations. (Default). 

=1, Calculations performed using 51 stations. (Recommended) 

ICOSZ - Interpolated Spanwise Sections Distribution Indicator 
=0, Equally Spaced. 

=1, Half-Cosine Distribution. 

IHEL - Helicoidal Surface Indicator: 

=0, Betap Surface used. (Default). 

=1, Helicoidal Surface used. 

PCTCL - (For short aft blades). Percentage of nominal span removed. 

Definition of Near Field Calculation 

INF - Near Field Calculation Option Indicator. 

=0, Far Field Calculation only. (Default). 

=1, Near Field Calculation under Far Field Conditions. (No fuselage or reflecting 
surfaces present). 

INFS - Near Field Calculation Spherical Limit Option. 

=0, Limits at Centre of Rotation, (ie; Far Field) 

=1, limits at Hub Radius. (Default) 

AKINF - Empirical Constant in Near Field Calculation; Default Value = 0.35 

Definition of Rotor 1 Wakes/Vortices 

IRW - Rotor Wake (i.e.. Rotor 2 Unsteady Loading Noise) Calculation Indicator. 

=0, No Wake Calculation performed. 

=1, Wake Calculation and Unsteady Loading Noise Calculation performed. 
(Default). 

IWAKE - Wake Model Indicator. 

=1, Linear Rational Function. 

=2, Kemp/Sears. (Default). 

=3, Mugridge/Morfey. 

ISHAPE - Wake Profile Shape Indicator. 

=1, Hyperbolic Secant. 

=2, Gaussian. (Default). 
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ITPVTX 


IHBVTX(*) 


IPRNTW 


NWHM 

SCD 

BETAW 

NSTEP 

VREF 

wnv 

Cl 

SBN(l) 

SBN(2) 

TAU 

tvti 


Tip Vortex Option Indicator. 

=0, No Tip Vortex Calculation. (Default) 

=1, Include Tip Vortex. 

Hub Vortex Option Indicator. 

=0, No Hub Vortex Calculation. (Default) 

=1, Include Hub Vortex. 

Wake/Vortex Velocity Profile Print Option Indicator. 

=0, No Printout (Default). 

=1, Print Profiles. 

Number of Rotor Wake Gust Harmonics Computed. (Default=Max.=21). 

NZ Values of Drag Coefficient 

Parameter to Account for Wake Flow Angle Variation from Freestream. 
(Default=0.0). 

No. of Points in Velocity Profile across one Blade Passage. (Max.=l50) 

Reference Velocity for Normalization of Gust Upwash Velocity Harmonic 
Coefficients, ft/sec. (Default=10.0) 

Inviscid Velocity Gradient, Normalized by Wheel Speed. (Default=0.0). 

Circulation Index for Tip Vortex Strength; (Default=1.0) 

Tangential Distance Normalized by Tip Blade Spacing for the Tip Vortex Center. 

Tangential Distance Normalized by Hub Blade Spacing for the Hub Vortex 
Center. 

Tip Clearance Normalized by Rotor Tip Chord. (Default=10.0) 

Tip Vortex Trajectory Index. (Default=2.0) 


NOTE: 

♦The Hub Vortex Option is present in the Code but is not calibrated. 


10 



Namelist ‘INSTAL’ 

(This File is only required if INST = 1) 

NPCLCD - Number of angles (max=5 1 ) at which Cl/Cd values are given (same for front and rear 
rotors) (integer input) 

NDES - Number of radial spanwise locations at which acoustic calculations are to be 
performed (Max=51) (integer input) 

NRADP - Number of radial spanwise locations at which pylon geometry is specified (max=5 1 ) 
(integer input) 

THETDP - Angular location of pylon: forward looking aft phi (degrees) 

CDP - Drag coefficient of pylon. 

IPYLON - IPYLON=0 will cause pylon effect to be not included in acoustic calculations 
(integer input) 

ALD - Angles at which front rotor Cl/Cd values are given: NPCLCD values are used 
(degrees): ALD(1).LT.ALD(2) LT.ALD(NPCLCD) 

CL - Front rotor lift coefficient values at angles ALD: NPCLCD values are used. 

CD - Front rotor drag coefficient values at angles ALD: NPCLCD values are used. 

ALDA - Angles at which rear rotor Cl/Cd values are given: NPCLCD values are used 

(degrees): ALD A( 1 ) .LT. ALD A(2) LT.ALDA(NPCLCD) 

CLA - Rear rotor lift coefficient values at angles ALDA: NPCLCD values are used 

CDA - Rear rotor drag coefficient values at angles ALDA: NPCLCD values are used 

RADP - NRADP radial spanwise stations at which pylon geometry is specified, normalized 
by R1 tip radius: RADP(1).LT.RADP(2) LT.RADP(NRADP) 

XLEP - Axial location of pylon leading edge relative to R1 PCA at RADP radii, normalized 
by R1 tip diameter. +ve upstream; NRADP values used 

CAMPD - Camber angles of pylon at RADP radii: NRADP values used (degrees) 

XTEP - Axial location of pylon trailing edge relative to R1 PCA at RADP radii, normalized 
by R1 tip diameter. +ve upstream; NRADP values used 

PSPANV - Percent span locations at which acoustic calculations are to be performed: NDES 
values are used: 0.00.LT.PSPANV.LT. 100.00: also PSPANV(1).LT.PSPANV(2) 
LT.PSPANV (NDES) 

THFDG - (representative thickness)/(max thickness) (used for thickness noise calculations) 

CLACC - Aircraft wing lift coefficient (0.0 = no wing present) 
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ARAC 

AATC 

WSS 

NPGC 

XWL 

YWL 

ZWL 

ARAD 

NHMAX 

NMODES 

NROT 

FLEN 

FDIA 

FRAWC 

NPHLAG 

ARBDU 


Aircraft wing aspect ratio 

Angle of attack of propellor re: freestream 

Aircraft wing semi-span (ft.) 

No. of pts. used to eval. Fourier coeff.’s (recommend 32) 

The following 3 inputs represent the X,Y ,Z location of the Rotor 1 disk center relative 
to the wing lifting line. 

X dist (ft) - positive downstream (aft) of lifting line 

Y dist (ft.) - positive horiz. along lifting line to a/c cl 

Z dist (ft) - positive vertically above lifting line 

Acoustic prediction arc radius (ft.) (overwritten by DIST if ISOLAT=l) 

Number of bpf harmonics for acoustic predictions (Max=5) 

Number of modes used in acoustic predictions (Max=10) 

Flag for front propeller rotation direction (FLA) -1 = clockwise 1 = counter 
clockwise 

Fuselage length (ft.) 

Max. fuselage dia. (ft.) (0.00 = no fuselage) 

Fraction of fuselage length ahead of wing lifting line 

Flag for phase lag in force response calculations 

0 = no phase lag 

1 = phase lag included 

Nondimensional once per rev axial velocity distortion (u/U) : negative on top (PHI= 
90 deg.) and positive on bottom (PHI = 270 deg.) for positive ARBDU 


The Installation-related input is printed in the installation output file, together with the results 
of preliminary aerodynamic calculations and the acoustic predictions as follows: 


* Preliminary Aerodynamic Output 


XV(1) - 

XV(2) - 

TOTTHF - 
TOTTHR - 
TOTTH - 
PWRF - 
PWRA - 
PTOT 


computed 70 % span pitch (front) after power match 

computed 70 % span pitch ( rear) after power match 

total thrust (front) 

total thrust ( rear) 

total thrust 

power (front) 

power ( rear) 

total power 
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♦Installation-Related Acoustic Predictions: 


I 

THTDIV 

PDBSL 

PHDV 

DIFFDB 


harmonic number index 

(THAD) input angles for acoustic predictions 

reference (clean) dB levels to which installation effects delta dB’s are applied. 
(PHI) azimuthal angles for acoustic predictions 

installation effects delta dB as a function of THTDIV and PHID, relative to reference 
dB levels at each THTDIV. 


SUBROUTINE GAMCAL (Page 54) 


Specific heat ratio calculation for air at low pressures from Keenan and Kaye gas tables. 

SUBROUTINE PITCH1 (Page 54) 

Calculate blade coordinate changes due to pitch angle change relative to ’design’ setting — also 
calculate Face Alignment (FA) and Mid-Chord Alignment (MCA) (Rotor 1) 

SUBROUTINE LSPFIT (Page 114) 

This routine is a least-squares parabolic curve fit program that will integrate or interpolate using 
a parabola which passes through points i and i+1 and misses points i-1 and i+2 (if they both exist) 
such that the square of the deviation is a minimum. 

Note that i is generally selected such that; 

x(i).le.xc.lt.x(i+l) 

the equation for the parabola is: 

y-y(i) = b*(x-x(i)) + c*(x-x(i))**2 

Outside of the x,y-data range, linear extrapolation of the parabola end point slope is employed 
if nxtrp = 0. 

SUBROUTINE PITCH2 (Page 56) 

Calculate blade coordinate changes due to pitch angle change relative to ’design’ setting — also 
calculate Face Alignment (FA) and Mid-Chord Alignment (MCA) (Rotor 2) 

SUBROUTINE RWVM (Page 65) 

This Subroutine computes the Rotor 1 Wake/Vortex characteristics and evolution with distance 
downstream. The following Subroutines are called from RWVM: 

VORTX1 - Computes strength and radius of tip (and hub) vortices (Page 72) 

WAKE1 - Calculates Wake Centerline Defect and Semi-Wake width (Page 70) 

WAKE2 - Calculates Tangential Wake Profile (Page 71) 
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V0RTX2 - Computes velocity field induced by vortices at all radii (Page 73) 

HRMONIC - Computes harmonic content of Rotor 1 wake/vortex flow at 1/4 chord point of 
Rotor 2 (Page 77) 

These routines are described in NASA CR-1815135, “An Investigation of Counterrotating Tip 
Vortex Interaction,” by R.K. Majjigi, K. Uenishi, and P.R. Gliebe (Reference 7). 

SUBROUTINE SKCAL (Page 64) 

This subroutine calculates a compressible form of the Sears Function (Unsteady Lift Response). 
The output is used in the Rotor/Rotor interaction noise calculation. 

The following Subroutine is called from SKCAL: 

FRESNL - Computes Fresnel Integral Functions C(Z) and S(Z) (Page 63) 

SUBROUTINE BJFN (Page 63) 

This Subroutine computes J-sub-n(x) Bessel Functions. 

SUBROUTINE AKVCAL (Page 59) 

This Subroutine calculates the Thickness or Volume Displacement noise source chordwise 
non-compactness Factor k-sub-v — for the following analytical distributions: 

Point-Impulsive 

Rectangular 

Parabolic 

Triangular 

Generalized Trapezoidal 
Generalized Parabolic 

SUBROUTINE AKLCAL (Page 61) 

This Subroutine calculates the Steady Loading noise source chordwise non-compactness 
loading factor k-sub-1 for the following analytical distributions: 

Point-Impulsive 

Rectangular 

Linear 

Generalized Parabolic 
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SUBROUTINE TOBSPL (Page 57) 

This Subroutine computes 1/3-octave levels given the tone harmonic spectrum. 

SUBROUTINE ATMOSP (Page 62) 

This Subroutine applies atmospheric air attenuation corrections for Standard Day (77 deg. F and 
70 pet. rel. hum.) conditions from the Shields and Bass correlation to lossless one-third octave 
spectra. 

SUBROUTINE SPLOUT (Page 58) 

This Subroutine Calculates dBA and dBD weighted levels, and prints out 1/3-octave SPL 
tables. 


SUBROUTINE CRPIE (Page 80) 


This Subroutine calculates ‘installation effects’ on Counter Rotation Propfan Noise. 

It requires access to the ‘IMSL’ (Edition 10) Library programs: ‘CERFE, FFTRF, NEQNF, 
U4INF, UMINF 


In addition, the following Subroutines are called from CRPIE or are User-Supplied to the IMSL 
routines UMINF and NEQNF: 

FCNP - Used to find 70% span pitch angles 

FCN - Used to find angles of attack on the front and rear Rotors 

CLCDEV - Calculates Q and Cd on front and rear Rotors given alphas 
SRPIE - Calculates tone noise deltas from installation effects 


SUBROUTINE FCNP (Page 94) 

This Subroutine is ‘User-supplied’ to IMSL solver UMINF. It is used to obtain 70% span pitch 
angles for the front and rear rotors such that the calculated power absorbed by each rotor 
approximately equals that input. 

The following Subroutines are called from FCNP or are User-Supplied to the IMSL routine 
NEQNF: 

CLCDEV - Calculates Cl and Cd on front and rear Rotors given alphas 
FCN - Used to find angles of attack on the front and rear Rotors 

SUBROUTINE FCN (Page 92) 

This Subroutine is ‘User-supplied’ to IMSL solver NEQNF. It is used to solve for the 
steady-state angle of attack (alpha) at each radial station for the front and rear rotors. The following 
routine is called from FCN: 
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CLCDEV - Calculates Cl and Cd on front and rear Rotors given alphas 


SUBROUTINE CLCDEV (Page 93) 

This Subroutine evaluates the Lift and Drag Coefficients on the front and rear rotors given the 
corresponding angles of attack. The following Subroutine is called from CLCDEV: 

LSPFIT - General interpolation and integration routine 


SUBROUTINE SRPIE (Page 96) 


This Subroutine calculates the effects of the installation environment on the BPF plus harmonics 
thereof of the forward and aft rotors. The effects are calculated separately for each rotor. The 
following Subroutines are called from SRPIE: 


BJFNA 

FCOEFF 

PHLAG 

OUTPUT 


Calculates Bessel (J) Functions 
Calculates installation-related velocity distortions 
Calculates phase lag of force response 
Outputs results of calculation 


SUBROUTINE BJFNA (Page 117) 

This Subroutine calculates Bessel Functions of the first kind for non-negative integer order and 
real argument, using the recursion method. 

SUBROUTINE FCOEFF (Page 107) 

This Subroutine evaluates u,v distortions at the plane of a propeller disk due to angle of attack 
operation, wing+trailing vortex system, the presence of a mounting pylon upstream, and the 
presence of a fuselage. It also carries out a Fourier Analysis of these distortions in the ’phi’ direction. 
In addition to the IMSL routine FFTRF, the following Subroutine is called from FCOEFF: 

FUSEL - Calculates fuselage induced flowfield 

SUBROUTINE FUSEL (Page 111) .. ..... 

This Subroutine computes the parameters needed to evaluate the fuselage induced flow, based 
on a Rankine Solid model. The IMSL routine NEQNF with the following ‘User-Supplied’ 
Subroutine is called from FUSEL: 

FCNFU - Function whose zero is to be found for fuselage flow 


SUBROUTINE FCNFU (Page 112) 

This Subroutine is ’User-supplied’ to the IMSL routine NEQNF. It computes the function 
whose zero is to be found to evaluate the fuselage induced flow based on a Rankine Solid model. 
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SUBROUTINE PHLAG (Page 113) 

This Subroutine computes two dimensional, linearised, compressible, flat plate convected gust 
phase lag formulae based on Amiet’s low frequency theory where applicable and on a Fresnel 
Integral based formula otherwise. See Goldstein: ‘Aeroacoustics’: pages 138-139, eqs.(3.70) and 
(3.71) for these formulae. In addition to the IMSL routine CERFE, the following Subroutine is 
called from PHLAG: 

BJFNA - Calculates Bessel (J) Functions (Page 117) 

SUBROUTINE OUTPUT (Page 106) 

This Subroutine writes out the Final Results of the predicted effects of installation environment. 
These Results are represented acoustically by Delta dBs and are presented in table form for each 
harmonic of BPF as Delta dB vs. angle Theta and angle Phi. Results for the front rotor are written 
first, followed by the results for the aft rotor. The Final set of information combines the installation 
effect Deltas calculated in CRPIE with the BPF harmonic tones predicted in the CRPFAN main 
program. 
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CRPFAN Source Code 



COUNTER-ROTATING PROPFAN NOISE PREDICTION PROGRAM 
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Rotational Speed, cpm. 
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CHORD(+) — Blade Chord normalised by Tip Diameter, b/Dt . 
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- Doppler Shift Indicator: 

=0, No Doppler Shift. 

=1 , With Doppler Shift. 
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Rotor Wake ( ie Rotor 2 Unsteady Loading Noise) 
Calculation Indicator. 

=0, No Wake Calculation performed. 

=1, Wake Calculation and Unsteady Loading Noise 
Calculation performed. (Default). 


DUB 1 : [ CHARLOTTE . NASA . DELIVER . CRP ] CRPFAN . FOR ; 5 9-OCT-1989 14:28 Page 
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SBN(2) — Tangential Distance Normalized by Hub Blade Spacing 

for the Hub Vortex Center. 
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THETDP - angular location of pylon : forward looking 

aft phi < degrees ) 
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the following 3 inputs represent the X , Y , Z location of the 
Rotor 1 disk center relative to the wing lifting line. 
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PHIT = PHIS+PHIO+PHIC 
NPH a PHIT/TPI 

PHITl a PHIT - TP I * FLOAT < NPH ) 
CSPHI1 = COS ( PHITl ) 

SPHI1 a SIN < PHITl ) 
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CALL LSPFIT(Z11 ,CHDU , NZ1 1 , RSEC1 , CHORDl , NRSEC 
CALL LSPFIT ( Zll , BP11 , NZ 1 1 f RSE C 1 , BETAP 1 , NRSEC 
CALL LSPFIT ( Zll , TM0C11 , NZ1 1 , RSEC1 ( TMOCl , NRSEC 
CALL LSPFIT (Zll , ZMCH , NZ1 1 , RSEC1 , ZMCl , NRSEC 
CALL LSPFIT ( Zll ( YMCll , NZ 1 1 , RS EC1 , YMC1 , NRSEC 
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1/3-OCTAVE LEVEL CALCULATION FROM TONE HARMONIC SPECTRUM 
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THIS SUBROUTINE CALCULATES THE THICKNESS 
DISPLACEMENT NOISE SOURCE CHORDWISE NON- 
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DETERMINE TOTAL PERTURBATION VELOCITY 


PRINT NUMERICAL VELOCITY PROFILES 
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300 CONTINUE 

GO TO (310,320), ISHAPE 

310 PP = 2 , / ( E X P ( 1 - 3169579* DT ) +• E X P { -1 . 3169579* DT ) ) 
GO TO 400 

320 PP = EXP ( — P I * DT * DT/4 . ) 
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'■O 

in 




* 


* 



.. 




« 


* 

Eh 






« 


* 

X 


m 




« 


* 

W 


w 




* 

U- 

* 

V) 


id 




* 

O 

* 

to 


U 




'* 


* 


t-4 

to 




41 

V) 

* 

— • 

ffl 

*> 




41 

D 

* 


» 

— «. 




41 

M 

41 

in 

O 

04 




41 

Q 

4t 

w 

o 

■v 




41 

< 

41 

Eh 

Eh 

O 




41 

04 

41 

w 

3: 

> 




* 


41 

X 


a 

H 


• — • 

4t 

a 

41 

Eh 

Q 

CO 

Eh 


Eh 

« 

z 

41 

CO 

04 

> 

> 


H 

41 

-< 

* 

CO 

O 

CO 

Eh 


x 

41 


<fc 

* 

X 

» 



U 

41 

X 

« 

— • 

u 


M 


CO 

41 

Eh 

* 


«. 

CM 

U 


* 

41 

O 

* 

in 

06 

-w» 



Eh 

41 

z 

* 

w 

X 

o 



< 

41 

U 

* 

3 

a. 

O 

i—4 


Eh 

« 

M 

41 

X 

j 

CO 

in 


W 

4t 

Eh 

41 

u 

< 


- — - 


X 

« 

Cfl 

41 

CO 

* 


06 



41 


« 

to 

D 

04 

ta 


VI 

41 

X 

41 

* 

< 

v 

< 


•» 

4l 

u 

41 

— . 

H 

o 



Eh 

* 

Eh 

41 

i-H 


U 

Eh 


< 

* 

as 

4i 

in 

X 

06 

O 


2 

* 

o 

41 


Eh 

M 

o 


— CD 

— > * 

~ 4 - 

o J > ~ 

< U < H M 

Eh — J m a, 

W H U Q * 

m oi ~ o ’ 

a &h < fN 

w 06 w os * « 

OS* O* o H 

* J o * U Q 

Eh << O > ffi O 


X X 

u u 

Eh Eh 

> > - 

\ V <N 


H r-i 

< - 

^ Eh — O 

U Eh 2 < X 

S £ < Eh O 

Eh U Eh td M 

VI V) < 03 W 


- U 'fl 'f H E 
\\«o h a* o 
< *< tn o * U 


Eh H 

Eh o o 

<0 U > 


— < 
O E 
{a. h dj 

M VI V) 


SSiflooPJoflpJEOH J 

U U • • • H ' O N O V) 4 Pi 

V)V)OOr\lUO<UO>UbH 


II II II If It If II fl II fi ff II I! II II II II 


•H — — — 04 

a Eh E-< O ~ O — -h 

H O « Q H SEho OV)Uh>h^O 

UE<2Bio«UWU >Qtf-Qw4^ 
SS:EHfl«Oos;2Qfr;FHQOMOV)t4{« 
EtUu4aECdE0HOWi<U0>U£t < 0 

V)lOCQ.<U!SCO><UO>tOCOlOV}V)Vl(J 


oooaooZ 

< 

o o o o o o * 


~ o -- o ~ r-4 06 

Q 06 — a — _J M 

OwowJrt — 
4jUO>UEufc 
to ui w ui w n h 


uu uuuuu u 


74 


SIGR = SSIGR(NSTR) 


DUB 1 : (CHARLOTTE. NASA. DELIVER. CRP ] CRPFAN . FOR ; 5 9 -OCT -1989 14:28 Page 57 
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2QRT = SDIST*C0S{ ALP'HR ) * R * SIG 1 *2 - * P I/B 1 
BRRT * TVTI *ZQRT/< 16 . 0 *ZQRT+i . 0 ) 

BR = BRRT * B 1/ ( 2 . *PI ) 

BR = BR + TAU * S IG1T 


DUB 1 : (CHARLOTTE .NASA. DELIVER. CRP ]CRPFAN. FOR; 
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FREE VORTEX NORMAL VELOCITY CALCULATION 
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INVL , LEND 4 NXTRP/0 , . FALSE . , . FALSE . 
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. EQ .1.0, CHK FOR INTEGRATION AND DOUBLE PT BEFOR INCREMENTING INVL 


SEARCH FOR APPROPRIATE INTERVAL 
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EXTRAPOLATION, ADJUST B, YI AND XD 


DUB 1 : ( CHARLOTTE . NASA . DELIVER . CRP ] CRPFAN . FOR; 
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360 I C = IC+1 

I F ( NXC- 1 C ) 380,370,370 
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t F ( Y . GT .8.0) GO TO 60 
Q=FM+5.0*(Y**0 .33333) 
GO TO 70 
CONTINUE 


DUB1 : [CHARLOTTE. NASA. DELIVER. CRP ] CRPFAN . FOR ; 5 9-QCT-1989 14:28 Page 101 
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Sample CRPFAN Input File: 


Rotor 1 Geometry 

+ 

Operating Conditions 




DUB 1 ; { CHARLOTTE . NASA . DELIVER . CRP ] F7_CRP_VTX . DAT ; 3 5-OCT-196 9 15:39 page 
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Sample CRPFAN Input File: 

Rotor 2 Geometry 
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Sample CRPFAN Input File: 


Installation Input 



DUBl : (CHARLOTTE .NASA. DELIVER. CRP]CRP_INST. INP; 1 5-OCT-1989 15:39 Page 
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Sample CRPFAN Output File: 


“Isolated” Prediction 




PROPELLER NOISE PREDICTION PROGRAM 
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TIP SPEED HACK NO. * 0.7114 
FLIGHT KACH NO. » 0.1998 
HELICAL MACH NO. = 0.7389 
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AIRCRAFT WING PARAMETERS: LIFT COEFFICIENT = 0-000 

ASPECT RATIO = 10-000 

SEMI-SPAN (ft.) » 6.800 
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SOME OVERALL CALCULATED STEADY STATE QUANTITIES 
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CALLS IMSL ROUTINES TO SOLVE C/R PROPELLER PERFORMANCE 
CL & CD VS ALPHA FORMULATION 
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AL1 , AL2 , CLA1 , CLA2 , CDAl , CDA2 
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THIS ROUTINE IS A LEAST-SQUARES PARABOLIC CURVE FIT PROGRAM THAT 
WILL INTEGRATE OR INTERPOLATE USING A PARABOLA WHICH PASSES THROUGH 
POINTS I AND 1+1 AND MISSES POINTS 1-1 AND 1+2 (IF THEY BOTH 
EXIST) SUCH THAT THE SQUARE OF THE DEVIATION IS A MINIMUM. 
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100 INVL = MAX ( 0 , MIN ( INVL , N+l ) ) 
WITHIN aa .FALSE. 
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LEAST SQUARE PARABOLA FIT FOR INTERVAL INVL 


DUB1 : (CHARLOTTE. NASA. DELIVER. CRP]CRPCLD. FOR; 2 5-OCT-1989 13:52 
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IllilkUil li ilti li III 11 liliiiiliil 



Sample CRPCLD Input File: 

Rotor 1 Geometry 
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Sample CRPCLD Input File: 

Rotor 2 Geometry 


/ CpQ' 
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Sample CRPCLD Input File: 
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Sample CRPCLD Output File: 


Printout 
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